Primordial black hole constraints on non-gaussian inflation models 
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We determine the abundance of primordial black holes (PBHs) formed in the context of non- 
gaussian models with primordial density perturbations. We consider models with a renormalized 
probability distribution function parametrized by the number, v, of degrees of freedom. We show 
that if 1/ is not too large then the PBH abundance will be altered by several orders of magnitude with 
respect to the standard gaussian result obtained in the v —> oo limit. We also study the dependence 
of the spectral index constraints on the nature of the cosmological perturbations for a power-law 
primordial power spectrum. 
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I. INTRODUCTION 

Contrary to early beliefs some level of non-gaussianity 
is expected to be present in all inflationary models. Fur- 
thermore, the non-gaussian contribution will in general 
be scale-dependent (see for example jl| and references 
therein). Consequently, a strong non-gaussian compo- 
nent on small-scales may turn out to be consistent with 
current observations which only allow for small devia- 
tions from gaussianity on large cosmological scales j^. 
Topological defect models also produce a scale-dependent 
non-gaussianity in a natural way However, in these 
models the fluctuations are not primordial but generated 
actively on increasing larger scales by defect network evo- 
lution. 

Black holes of various masses are expected to form in 
the early universe as a result of the collapse of large den- 
sity perturbations 0, . The abundance of PBHs pro- 
vides useful constraints on the primordial power spec- 
trum over a wide range of small scales (see for exam- 
ple d, S S IS Q and references there in) . Although, 
the usual calculation assumes gaussian perturbations (see 
however ref. |^ for a discussion of the role of non- 
gaussian fluctuations in PBH formation), standard con- 
straints will need to be modified when considering non- 
gaussian models. 

The cosmological constraints on the abundance of pri- 
mordial black holes come either from their contribution 
to the matter density at the present time or from the cos- 
mological implications which result from their Hawking 
evaporation PBHs with mass M < 10^^ will have 

evaporated by the present epoch but if they have a mass 
M > 10^ they will still be present at the nucleossynthe- 
sis epoch. Assuming a standard cosmological scenario, 
the fraction of the universe energy density in PHBs with 
mass greater than M at the time they form is constrained 
by observations to be < 10~^° on any interesting mass 
scale with more accurate constraints depending on their 
mass (see for example ref. 0|). 

We will generalize to the non-gaussian case the stan- 
dard calculation of primordial black hole abundances. 



This uses the Press-Schechter approach [I^ in order to 
estimate the fraction of the energy density in PBHs. An 
alternative calculation using peaks theory was performed 
in 9] using the criterion for black hole formation derived 
by Shibata and Sasaki which refers to the the met- 
ric perturbation rather than the density field. However, 
given the large uncertainties involved the results obtained 
by the two methods were shown to be consistent. 



II. THE POWER SPECTRUM 

On comoving hypersurfaces there is a simple relation 
between the density perturbation. A, and the curvature 
perturbation, TZc (e.g. 0) 
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Here w — p/p is the equation of state, H = d/a is the 
Hubble parameter, fc is a comoving wavenumber and c is 
the speed of light. The power spectra are related by 



and consequently at horizon crossing 
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If we assume a power-law primordial power spectrum 
with Vtz^ (k) — A-ji^ (k/ko)"'^ and choose a gaussian win- 
dow function, 



W{kR) = exp 



k^R^ 



(4) 



then the variance of the primordial density field 
smoothed on a scale R at horizon crossing. 
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is given by 
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for n > 1. We note that for n < 1 the variance of a den- 
sity field with a power-law primordial power-spectrum 
diverges due to contributions near fc = 0. However, a 
small k (large wavelength) cut-off to the power spec- 
trum is expected on physical grounds. Spergel et al. 
have shown, using the WMAPext-)-2dFGRS dataset, that 
Atz, = (0.8 ±0.1) X 2.95 x 10"^ for fco = 0.05Mpc"\ We 
shall take w = 1/S appropriate for the radiation domi- 
nated era. 



III. THE PBH ABUNDANCE 

We shall use the Press-Schechter approximation (0|) 
to compute the fraction of the energy of the Universe, 
/(> M), associated with PBHs with masses larger than 
M at the time they form. The Press-Schechter approx- 
imation was originally proposed in the context of initial 
gaussian density perturbations and was much later gen- 
eralized to accommodate non-gaussian initial conditions 

m 

The mass fraction is assumed to be proportional to 
the fraction of space in which the linear density contrast, 
smoothed on the scale M, exceeds a given threshold Ath: 
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Here Pfl(A) is the probability distribution function 
(PDF) of the linear density field A smoothed on a scale 
R and ^/ is assumed to be a constant which can be cal- 
culated by requiring that /(> 0) = 1, thus taking into 
account the accretion of material initially present in un- 
derdense regions (note that Af = 2 in the case of gaussian 
initial conditions). 

This generalization has successfully reproduced the 
results obtained from A'^-body simulations with non- 
gaussian initial conditions [l^ but it has been shown in 
ref. 10 that it does not adequately solve the cloud-in- 
cloud problem (see also ref. ^^)- Although this means 
that in many models ^/ is not expected to be exactly 
constant, it will in general be of order unity. Conse- 
quently this small uncertainty will have a negligible im- 
pact on our final results and we shall not consider it 
further in the following analysis (we will use Af = 1 
throughout the paper). 

The precise value of the threshold value Ath relevant 
for the calculation of the PHB abundance is uncertain 
but we shall use the standard value, Ath — 1/3, for a 
radiation dominated universe . Although there should 
also be a finite upper limit to the integration in Eq. Q 
in practice the upper cut-off is unimportant since Pr(A) 
is usually a rapidly decreasing function of A. 



In order to implement our calculation one needs to re- 
late the mass M to the comoving smoothing scale R. An 
overdense region in order to collapse must be larger than 
the Jeans length at maximum expansion which, in the 
radiation dominated era, is 3"^^^ times the horizon size. 
On the other hand, it cannot be larger than the horizon 
size or otherwise it would form a closed universe separate 
from our own _21J. For simplicity, here we will assume 
that PBHs will form with roughly the horizon mass. This 
assumption is only approximately true (in general the 
mass of primordial black holes is expected to depend on 
the amplitude, size and shape of the perturbations). It 
will nevertheless have a small impact on our main results. 
Hence, when the comoving scale R = c{Ha)~^ enters the 
horizon the horizon mass is 
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and primordial black holes with a mass M = Mh may 
form at that time. In the radiation dominated era p oc 
5* ^^^a^^ where is the number of relativistic degrees of 
freedom (we approximate the temperature and entropy 
degrees of freedom as equal) . This implies that 
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Taking into account that fl^^oh^ = 2.5 x 10 ^ it is 
straightforward to show that 
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The number of relativistic degrees of freedom, g^,, in the 
early universe is expected to be of order 100. However, 
the smoothing scale R is only weakly dependent on g^, 
(note that Eq. H10|l corrects a minor error in ref. |^). 



IV. THE X MODEL 

We consider an initial density field with a chi-squared 
probability distribution function (PDF) with z/ degrees of 
freedom, the PDF having been shifted so that its mean 
is zero (such a PDF becomes gaussian when v oo). 
Hence, 



/(> M) = Q{u/2, Ath/<7A X vAV2 + 1^/2) 



(11) 



where Q{a, x) with a > is the incomplete gamma func- 
tion defined by 
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FIG. 1: PBH abundance as a function of black hole mass for 
power-law power spectra with n = 1.2, n = 1.3 and n = 1.5 
(sets of curves from left to right respectively) in a model 
with n — oo, n = 100, n = 10 and n = 1 degrees of freedom 
(solid, dashed, dot-dashed and dotted lines respectively). 

and r(a) = r(a,0). If — s- oo then 

We further assume that the shape of the PDF is scale 
independent, that is Pr{A) = P{R, A/aA{R))/(JA{R) is 
always the same function (or equivalently v does not de- 
pend on R) . Ahhough this is not expected to be precisely 
true in realistic models, it may be a good approximation 
over the range of scales relevant for the determination of 
PBH abundances. 



V. RESULTS AND DISCUSSION 

In Fig. n] we plot the PBH abundance as a function 
of black hole mass assuming a power-law spectra with 
n = 1.2, n — 1.3 and n = 1.5 (sets of curves from left 
to right respectively). The gaussian model is represented 
by a solid line and the dashed, dot-dashed and dotted 
lines represent models with n — 100, n — 10 and 
n = 1 degrees of freedom respectively. We clearly see 
that even a small amount of non-gaussianity may have 
a large impact on the estimated PBH abundance. This 
shows that an improper use of the standard gaussian as- 
sumption may introduce errors on the estimated PBH 
abundance of many orders of magnitude. 

It is useful to determine the impact of non-gaussianity 
on the spectral index constraint, n. We clearly see from 
Fig. ^that for x"^ models with a finite numbers of degrees 
of freedom the constraints are strengthened since for the 
same value of the spectral index the PHB abundance is 
significantly increased. Still, this will at most shift the 
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FIG. 2: The value of a (defined as the ratio between the 
amplitudes of the perturbations in a standard gaussian model 
and in a x^ model with identical mass functions, /(> M), at 
a given mass, M) as a function of v and /(> M). Note that as 
expected a decreases with increasing v (approaching 1 when 
v — > oo). Also, for a fixed v a smaller /(> A4) implies a larger 
a. 

spectral index constraint down by about 0.1 (in the most 
extreme case with v = I). 

It is also interesting to ask how much larger the ampli- 
tude of the perturbations has to be in a standard gaus- 
sian model relative to a model in order for the mass 
functions /(> M) at a given mass to be the same. An 
answer to this question is easily obtained by solving the 
equations 

ierfc(x) = /(>M), 
Q{i^/2,xaV^+iy/2) = f{>M). (15) 

with respect to a for a given ly and /(> A/). The results 
are shown in Fig. |21 As expected we clearly see that 
if v is increased then a decreases (approaching 1 when 
1/ — > oo). Also, for a fixed a smaller /(> M) implies 
a larger a since the differences between gaussian and x^ 
probability distribution functions (for v ^ oo) become 
more dramatic as we move deeper into the distribution 
tail. 

In this study we do not compute the PBH for any spe- 
cific model of infiation. Still, we expect the class of 
models we investigate in this paper to be representative 
of a large family of non-gaussian models which may be 
realized in the context of specific infiationary scenarios. 
The simple assumptions made in this paper also have 
the advantage of allowing for a complete separation be- 
tween the effect of the power spectrum (which controls 
the amplitude of density perturbations on the scales rel- 
evant for the calculation of PBHs abundances) from the 
distribution of the phases of the various Fourier modes. 
Although, it is possible to consider stronger deviations 
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from a gaussian random distribution which would mod- POCTI/FP/FNU/50161/2003. 
ify the constraints even further, it might be difficult to 
construct realistic models where such deviations are in- 
corporated in a natural way. 
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